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investigations, incJuding studies of crustal dynamics, and earth 
and ocean tides, and for the general development of precision 
orbit determination. The network performed tegular tracking of 
several other retroretlector satellites including GEOS-1, GEOS-3, 
BE-C, and Starlette for refined determinations of station 
coordinates and earth's gravity field and tor studies of solid 
earth dynamics. 

A major program in laser upgrading was initiated during this 
period to improve ranging accuracy and data yield. This program 
included 1) an increase in pulse repetition rate from 8 ppm to 30 
ppm, 2) a reduction in laser pulse width from 6 nsec to 2-3 nsec, 
3) improvements in the photoreceiver and the electronics to 
improve daylight ranging, and 4) an analog pulse detection system 
to improve range noise and accuracy. 

Tests were begun at Mt. Hopkins on three candidate analog 
pulse detection systems to replace the waveform digitizers 
(WD2000) . The major impetus for the analog system is the speedup 
in data rate, which could not be accommodated with the digitizer 
and the increased difficulty in maintaining the complex digitizer 
equipment. The analog system will also offer potential for 
improved range noise and accuracy. A design was completed for 
the modifications to the photoreceiver for the narrow band filter 
and a more flexible arrangement for the PMT Mounting. A new PMT 
and PMT base were also selected. 

During the reporting period, work continued on the 
replacement for the bi-planar photodiode currently being used in 
the laser system start channel. The PIN diode should give a 
faster, more stable response (less pulse distortion) and less 
jitter in the system start circuit. 


PAGE 4 


In May, a visit was made to Bangalore, India, to discuss the 
relocation of the Mt. Hopkins laser system with officials of the 
Indian Space Research Organization ( ISRO) . There appears to be 
agreement on locating the laser at Naini Taij however, the 
response from ISRO has been very slow. An updated draft 
Memorandum of UnderstandirT was sent to ISRO from SAO. We still 
are awaiting their response. Tn November the NASA Administration 
discussed the Laser Program with ISRO while he was in Bangalore. 
Dr. Dhawan agreed to provide a response in early CY 1981. 
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2. OPERATING STATUS 

The SAO laser sites in Natal, Arequipa, and Orroral Valley 
continied routine operations throughout the six-month period. 
Together with the cooperating stations in Wettzell, Shanghai, 
Dionysos, Grasse, Kootwijk, Helwan, Metsahovi, and Dodaira (see 
Section 3.), the laser stations obtained a total of 56,262 
observations on 3,341 passes of GEOS 1, GEOS 3, BE-C, Starlette, 
and LAGEOS. These statistics are quick-look passes and points 
and represent only a percentage of the actual data taken. 
Monthly statistics of the passes and points, by station and by 
satellite, are given in Table 1; these data have been sent to 
the National Space Science Data Center at Goddard Space Plight 
Center (GSFC) . Table 2 shows the 1980 final data processed to 
date. 


The Mt. Hopkins site was closed for regular data taking 
during the reporting period. The station remains prepared to 
provide tracking support for brief periods, if requested, but 
current activities are concentrated on upgrading, documentation, 
and preparation for relocation. 

We continue to maintain the operations reporting procedures 
requested by NASA by providing statistics of tracking success, 
weather, and maintenance on a monthly basis. Table 3 gives the 
six-month summary of this information. 
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Table 1. 

Quick-look passes and points, 1 July through 31 December 1980 
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Laser operations summary, 1 July through 31 December 1980. 
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3. LASER OPERATIONS AT COOPERATING AGENCIES 

In addition to operating its own stationsr SAO supports the 
operation of overseas cooperating laser stations by furnishing 
orbital elements for predictions, screening quick-look data, and 
acting as a general U.S. interface. These cooperating agencies 
are located in Greece, Japan, Spain, The Netherlands, France, 
West Germany, Egypt, Poland, and Finland. During the past six 
months, SAO actively supported the Centre National d* Etudes 
Spatiales (CNES) Starlette program with routine laser tracking, 
providing CNES with orbital elements to sustain its tracking 
operation. 

The Delft University of Technology’s laser system in 
Kootwijk, The Netherlands, obtained large amounts of tracking 
data on retroref lector satellites throughout the past six months. 
For this effort, SAG supplies orbital elements for laser 
predictions, provides routine processing of quick-look data, 
helps coordinate operational activities, and routinely screens 
the data. The laser has been operating at an estimated accuracy 
of 15 cm, with additional upgrading being planned. 

The laser system operated under the cooperative laser 
tracking program, begun in 1977, with the Soviet Academy of 
Science, the Technical University of Prague, Czechoslovakia, the 
Helwan Observatory in Egypt, and SAO, continued to provide data 
for research programs early in the reporting period. For this 
progr2un, the Soviet Union and the Czechoslovakians provide and 
maintain the laser tracking system at Helwan, the Helwan 
Observatory furnishe j personnel to operate the system, and SAO 
supplies technical consultation, a station clock, and partial 
operating support from PL 480 funds (Smithsonian Institution 
Foreign Currency Program) . SAO routinely screens and validates 
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the data. Dutiny the laat yeat, the system was upgraded with a 
new short-pulse laser, improved software, and more rugged optics. 
The improved range accuracy is estimated to be about 25 cm. 

The laser site in Borowiec, Poland, was not operational 
during this reporting period. This station was also operated 
under partial support from the Smithsonain Institution Foreign 
Currency Progreim. 


The station in Dionysos, Greece, owned by the National 
Technical University (NTU) , was not operational as it is still 
undergoing an upgrading program for improved ranging capability. 
This system should be back in operation in July in preparation 
for the MERIT Campaign. 

SAG also has close working arrangements with the If'Stitut 
fur Angewandte Geodasie (IFAG) in the Federal Republic of 
Germany. The IFAG satellite-ranging system in Wettzell has been 
tracking LAGEOS and other laser retroref lector satellites with 
its short-pulse neodymium Yag mode-locked laser system which has 
consistently demonstrated range-noise performance of 2 to 4 cm. 
The station operated during January, February and early March. 

A new site in Metsahovi, Finland, began sending data during 
this reporting period. The tracking system uses a 25-nanosecond 
Q-switched laser pulse with an estimated accuracy of 1 meter. 
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4. SATELLITE OBSERVING CAMPAIGNS 

The laser tracking network continued its program of data 
acquisition, with particular emphasis on tracking support for 
LAGEOS for the preliminary MERIT Campaign. In addition, 
satellite observations were made to: 

A. Support the scientific and crbital 
maintenance requirements for LAGEOS and the Laser 
Earth Dynamics Project. 

B. Support the study of earth body and ocean 
tides, seasonal and other variations in the earth's 
gravity field, and the investigation of polar motion. 

C. Provide data for improving the accuracy of 
station coordinates and the gravity-field model, which 
are necessary for LAGEOS and other geophysics 
programs. 

D. Support the tracking campaign for Starlette 
in conjunction with CNES. 
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The practical and scientific importance of the monitoring of 
the rotation of the Earth has long been recognized. The 
International Latitude Service was established in 1895 to study 
polar motion (i.e., the motion of the axis of rotation with 
respect to an axis "fixed" in the Earth ) , and has been succeeded 
by the International Polar Motion Service, now based in Japan. 
The Bureau International de I'Heure was established in 1911 to 
provide an international time-scale based on the rotation of the 
Earth, and today it is responsible for the international scales 
of atomic time and of universal time, and hence for the basic 
data on the variations in the rate of rotation of the Earth. An 
accurate knowledge of universal time and polar motion is requ’ •'ed 
for practical purposes, such as geodetic surveying and prec ue 
navigation, but it is also of great scientific value since it 
provides information about the interior of the Earth and 
geophysical phenomena that cannot be obtained in other ways; 
this information could itself be of practical application in the 
studies of earthquakes and climatic variability. 

At present, the international services rely almost entirely 
on data provided by a network of cooperating observatories using 
optical instruments such as zenith telescopes and astrolabes. 
New techniques, such as laser ranging and radio interferometry, 
promise to provide an order of magnitude improvement in the 
precision of the data, but (apart from polar-motion data obtained 
from the doppler tracking of active satellites) there have been 
no international arrangements for regular observations, nor for 
the collection and analysis of the data, nor for the speedy and 
regular publication of results. 

MERIT is a program sponsored by the International 
Astronomical Union (lAU) and the International Union of Geodesy 
and Geophysics (lUGG) for a special period of international 
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collaboration in the Monitoring of Earth-Rotation and in the 
Intercomparison of the Techniques of observation and analysis. 
The objectives of the project are: (1) to foster the development 
of new techniques for the measurement of the variations in the 
rate and axis of rotation of the earth, (2) to obtain precise 
data on Enrth rotation in order to increase our understanding of 
the causes and effects of the variations, and (3) to make 
recommendations on the observational basis and organizational 
arrangements for future international services on Earth rotation. 
The program is focussed primarily on LAGEOS for the measurements 
of polar motion and earth rotation, with Starlette and GEOS-3 
being tracked at lower priority for other investigations in 
geophysics and orbital analysis. 

The full proposed program of activities includes (1) initial 
short periods of observation (during 1900) to test the techniques 
and improve the arrangements for international cooperation; (2) 
a later common period (at least a year 1983/84) of regular 
observation by all suitable techniques; and (3) periods of 
planning, data analysis, and review. The program relies on the 
support of observatories, data analysis centers, and 
international organizations to ensure useful scientific results 
and a sound basis for the future development of the international 
services on Earth rotation. 

SAG has participated in MERIT in two ways. First, as a data 
center SAG has provided predictions, status reports, and 
collected observations for the Satellite Tracking Laser Systems. 
Second, as an analysis center SAG has used the quick look data 
received each week, and calculated the position of the earth's 
pole and the Length of Day (LGD) . These data are supplied weekly 
to the Central Coordinating Center at the Bureau International de 
I'Heure (BIH) . At present SAG is the only new group to provide 
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these data routinely in an operational way. 

The accu.acy of the SAG pt osition and LOD is hard to 
assess since there is, at present, no suitable standard of 
comparison. We believe it to be an improvement over classical 
data . 


To prepare for the three-month Preliminary MERIT campaign, 
SAG began routinely determining polar motion from Quick Look Data 
from January 1980. We now have almost one year of data. As a 
byproduct, we have established the coordinates of many new Laser 
stations which participated. During this period alone, 20 laser 
stations provided data, giving one of the richest data sets ever 
obtained. These coordinates seem to be known with an accuracy of 
several decimeters, and continued refinement of these coordinates 
during the next few years may also result in measurement of 
relative motion of the observating stations, as well as the 
motion of the whole earth in space. 

From the success of the MERIT campaign and examination of 
this year's data from SAG, both BIH and the U.S. Naval 
Gbservatory are planning to use the SAG polar motion and earth 
rotation data in the officially adopted data. SAG intends to 
continue this Quick Service and may be a part of the new time and 
pole position service sanctioned by the International Unions. If 
so, this would be the first of the new space techniques that will 
have been incorporated into the classical astronomical services. 

Even aside from MERIT, the major emphasis of the SAG program 
has been and will continue to be on LAGEGS, with the focus on 
maximizing data yield and accuracy. 
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LAGEOS serves as a stable reference platform for 
observations in support of many investigations, including 1) the 
relative movement and deformation of the earth's tectonic plates, 
2) variations in the motion of the earth's polar axis and 
rotation, 3) improved geodetic and astronomical reference 
systems, 4) solid-earth and ocean tides, and 5) satellite orbital 
perturbations . 

SAG, GSFC, the University of Texas and several other 
organizations are using LAGEOS data to develop new analytical 
techniques for improved baseline measurements and earth-dynamics 
investigations. In particular, translocation techniques 
combining long-arc and overlapping observations from station 
pairs and other combinations of dynamic and geometric techniques 
offer promise of reaching the centimeter accuracies required to 
observe crustal motion. 

SAO lasers continued to track LAGEOS both for orbital 
maintenance and for geophysical and orbital analysis at SAO, 
GSFC, and other research organizations. With installation of the 
pulse choppers, the SAO lasers are now ranging to LAGEOS with an 
estimated systematic accuracy about 10 cm (see Section 8) . 
Returns are in the range of 1 to 10 photoelectrons, however, and 
shot-to-shot range noise is typically 30-50 cm owing primarily to 
photon quantization and sampling errors. To improve the quality 
of the data, SAO continues to examine all systematic influences 
at these low signal strengths and to tighten the screening 
procedures. During this reporting period, a program of laser 
upgrading has been undertaken to improve data accuracy and yield 
on LAGEOS (See Section II) . 
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A summary of the network's tracking activity on LAGEOS 
during the reporting period is included in Table 4. 

Table 4. 

Quick-look LAGEOS passes 
1 July through 31 December 1980 


Station July 

Arequipa 247 
Natal 60 
Orroral Valley 128 


435 


Aug . 

Sept . 

Oct 

352 

302 

199 

147 

85 

125 

161 

164 

148 

660 

5 51 

472 


Nov . 

Dec . 

Total 

211 

161 

1,472 

137 

82 

636 

78 

64 

743 

426 

307 

2,851 


Total 
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5. ENGINEERING 

The Engineering Group of the Experimental Geophysics 
Department provides the daily hardware and systems support 
necessary to maintain routine network operations. It is also 
responsible for the system modifications and improvements 
required for new programs. 

5.1 Lasgc F h9t ..Q.re££ive]:. 

The coatings of the primary and secondary mirrors in the 
Natal site photoreceiver had deteriorated significantly in the 
past years of operation because of ambient salt air exposure. 
The opportunity was taken to strip, recoat and install the 
mirrors from the he. Hopkins photoreceiver into the Natal system 
resulting in an improved signal strength of approximately 3. The 
mirrors of the Natal system were in turn refurbished and 
installed in the Mt. Hopkins system. At that time the 
Mt. Hopkins photoreceiver was also upgraded to eliminate 
mechanical instabilities in its field-stop and flip-mirror 
mechanisms that had been problem sources since its prototype 
days . 


As a part of the laser upgrading program (see Section 10) a 
redesign of the photoreceiver has been developed to incorporate a 
3 Angstrom filter and a PMT shutter. This new design will enable 
the instrument to accommodate different PM tubes of varying 
photocathode sizes. The components needed for the modification 
will be purchased in early CY 1981. As a part of this upgrading 
program tests are underway to choose a better PMT and PMT base to 
reduce pulse distortion and range abasing error. 
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A iiew cooling unit purchased for the Mt. Hopkins laser was 
integrated and field tested at that site, making that system more 
reliable in future operation. 

In order to imp.^ove the lifetime of optical components in 
the laser transmitter, new calorimeters were purchased to assist 
in system calibration. These units were tested and found to be 
unstable while also producing apparent low readings. The 
instability problem was cured by replacing some marginal 
components but the question of low readings awaits further study. 
We expect to test the calorimeters again in the near future with 
a ruby laser in the Center for Astrophysics Shock Tube Lab as 
soon as that system becomes available to us. 

Since the introduction of the pulse choppers, it has been 
noticed that the laser start pulse is subject to noise and 
pulse-to-pulse variation. The noise source was found to be 
Pockels Cell switching which is picked up by the coaxial lines of 
the detector. The pulse variation has been attributed to pick-up 
in the detector and its pulse forming circuitry. To address 
these problems, we plan to cover the Pockel's cells with RF 
shields for the cells, double shielde the coax lines, and replace 
the current planar photodiodes with PIN diode detectors. Several 
versions of PIN diode detectors were tested and evaluated. To 
utilize these detectors, fast pulse amplifiers were purchased to 
increase the low output of the diodes to the needed voltage 
levels. Acceptance testing of the amplifiers showed them to be 
out of specification and they were returned to the manufacturer 
for repair. 

Studies to increase the laser transmitter repetition rate to 
30 ppm have been completed. A vendor. Laser Industries, has been 
selected to furnish the components required for modifications to 
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the laser power supply. See Section 10 for more detail. 

Efforts have continued to reduce the jitter in the return 
circuit pulse processor electronics and to construct an analog 
detection system to replace the WD2000 waveform digitizer. The 
analog detector would enable us to simplify the system 
electronics and is a fundamental step in the speed up of the 
laser repetition rate (see Section 10). 

Prototypes of three different analog detectors were built 
and tested in Cambridge with promising results. The prototype 
systems include: a system based on constant fraction 
discriminators, a distributed delay line system with amplitudes 
summed and differenced to measure center of area, and a 
combination of matched gaussian filter/pulse differentiator and a 
cross-over dectector. Testing of these devices has now moved on 
to Mt. Hopkins. 


5.3 MinicQmputers 


Each of the SAG laser field stations has a minicomputer 
system that performs system calibrations, ranging operation, and 
other auxiliary functions. Magnetic tape (Line tape) is the 


primary storage medium for 
interaction with the computer 
(Miniterm) and a CRT display. 


programs and 
is through a 


data. Operator 
keyboard/pr inter 
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The minicomputers and associated peripheral devices in 
Natal, Arequipa, and Ori'oral Valley continue to function 
normally. The system at Mt . Hopkins remains operational. 

All field minicomputers remained operational during this 
reporting period. A set of Line Drives was received from Peru 
for repair but that system remained operational with a spare set 
of drives. Repair work on the faulty drives is underway. 

In addition to the field minicomputers, we maintain four 
similar units in Cambridge. System A is used to post-process 
field station Line tape data, which is then stored on 
industry-compatible magnetic tape for subsequent transfer as 
input records to the VAX 11/760 central computer. System B is 
used in the development and testing of software for the field 
systems. System C serves as the testing mainframe and as a 
back-up in case of failures in one of the other systems. The 
fourth system is utilized routinely by personnel in the 
Communications division. 

Because of continued failures experienced on system B, 
system C was rewired to replace system B in order to more fully 
investigate this problem. A major failure occurred on system A 
which caused the loss of many disc files. The problem was 
isolated and corrected, but work continues on file recovery. 

The communications minicomputer continues to assist 
communications personnel in the preparation, storage, and 
archiving of data messages. Since the implementation of the data 
links, laser data have been routinely transfered to and from the 
VAX with a high degree of reliability. This computer is equipped 
with Line-tape units, a keyboard/printer, a paper tape reader and 
perforator, and a Lexiscope display. The iinframe and all 
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peripherals have been operating routinely for the last six-month 
period . 

A nine-track tape drive for one of the Cambridge based 
minicomputers has been received. This device, once interfaced, 
will eliminate a large part of the formatting and data transfer 
problems that are now experienced with the seven-track drive. 

Studies continue on the feasibility of utilizing a 
minicomputer-radio link-minicomputer communications path as well 
as further automation of the the Communications Center. These 
studies are aimed at increasing current system flexibility and 
reliability as well as accommodating any further expansion. 

5.4 T imekeeping 

During the last half of 1980, timekeeping systems for the 
SAO tracking network have maintained epoch time traceab.'e to (JTC 
(U.S. Naval Observatory) with an accuracy of better than plus or 
minus 3 microseconds, except for Australia, which maintained 4-/-5 
microseconds from July 1st through July 22nd, and Egypt, which 
maintains time to plus or minus 50 microseconds. Each SAO 
tracking site is equipped with a broad-based timing system 
comprised of dual parallel tim'.i-g channels. Cesium oscillators, 
backed up by rubidium oscillators, offer a stable time base for 
each channel. Redundant time accumulators guard against time 
discontinuities, and redundant VLF/OMEGA receivers provide a 
reliable backup and frequency reference for the system. Without 
NTS epoch data, portable clock comparisons are required to 
provide the necessary epoch reference checks until another 
satellite-based time transfer system, such as GPS or Transit, can 
be implemented. Table 5 summarizes the network accuracy. 
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All laser sites have had portable clock comparisons during 
this period, with the exception of Egypt. Table 5 summarizes the 
results of these comparisons 

The U.G. Coast Guard loaned SAG a portable cesium clock for 
comparisons at South American tracking stations, to prepare for 
the MERIT tracking campaign. Both the Natal and the Arequipa 
comparisons made in September tracked to the last comparisons 
made in June, 1978. Traceability was better than +/-4 
microseconds for Natal. The accuracy of the portable clock set 
was better than 0.4 microseconds. An earlier attempt at a 
portable clock trip in July failed because the power 
supply/charger developed problems ano only a one way comparison 
with the Observatorio Nacional in Rio de Janeiro could be made. 
Repairs and tests were made by the manufacturer and at SAG before 
the successful trip was completed in September. 

A portable clock trip to Grroral Valley during September 
helped explain epoch time differences during the June thru August 
period. Cesium replacements at the NASA Minitrack site during 
June and August caused some unexplained differences of up to 28 
microseconds between Australian National Mapping and NASA timing 
data. SAG derives its epoch reference from the minitrack site. 

A portable clock comparison was made between the Mt . Hopkins 
timing system and one at nearby Fort Huachuca, Arizona. Accurate 
time was necessary at the tracking site because of electronic 
systems tests to be performed starting in December. 

New 50 ohm driver boards were installed at both South 
American sites. Timing system counter measurements have improved 
with the better noise rejection and crosstalk characteristics of 
the new boards. 
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Line 10 TV' tirring equipment was installed at the Natal site 
to provide a daily epoch check. These comparisons are made among 
other sites in Brazil and are coordinated by the Observatorio 
Nacional in Rio de Janeiro. 


SAC NETVkCKK TIMEKEEF 


ii 

) 


■■n 

> 


j 

cx> 

cn 







C 3 

,i -3 



fi 



J 

0 

o 

CO 



Eh 





V 4 . 








iJ 

Vj 

HJ 












a: 



0 

‘ J 

•d 

t: 



>: 



iJ 

t 3 


* > 






0 

U 


r-( 



u 



-J 

w 


+3 >1 

rH 






o 

; 

c* 

i .3 a 






t: 

C 

c: 

>Jt o 





o 

0 

Vi. 


o 

o 




03 


r-- 


03 

03 

o\ 


o 







t 


CO 

iJ 

V-. 

'j 



iJ 


C.' 


u- 


03 


o 

cll 


O c 

iJ 

03 

>: 

a 



03 



Ot 

CO 





to 

lu 

,)k« 

w 


t i 

o 


VO 



.ii 

"JJ 


(*) 

(0 


-s 

rH 

'T 

< 


rH 





f-H 

r H 

tM 









30 



r ^ 

rH 


rH 





c. 


n 


n 


c,-; 



o 

• 4 -> 

ij 


vJ 

i .3 

*v' 





X 


X 


■ ^ 

E" 


;"’t 




• ^ 


O 

F t 

CO 

< 





^rV 



;j 

If'} 

< 

< 


< 

< 

CO 




C/i 

t/i 


CO 

Ji 

» .. 


a 

tv? 


«/) 

. j 

E-* 

< 


’O 

Tm 

Eh H *j 





a 

rj 



f/j rH 

w 












fi 


c 




rC 

u 

rH 











rj 


H 

u 



3-1 £h 

H 



rH 



CT» 

Cj t j 

!— J 







O 



fj> o; 

V 







f 4 







fO 

G 



u 3 

rH 














,i.f 



I'l 

a. o 

u 













V", 

a 

■«: 


U 

u 3 

1 












O to 


o 


G 















5 »J 

'G 

*,y 


U 


f 











.y 

-13 


c; 

»*H 


G 













o 




G 















u 

fTJ 

l>-t 

' } 



■J 















M 

'O 

U .3 


) 


v'/ 











u 

.•«; 

Vi? 

0 ? 

v-^ 




•'O 












Ti 

iJ 

u 

;3 


u 

e-i 

c 











t: 

'D 



iij 


li) 

O v; 

o 











rH 

rH 

> 


G 


‘-M 

► H KH 

G 









ro 

rn. 

ro 

'J 

u 

u 

Cii 


if} 

Eh < 

w 

»r> If .' *x» 

o ij r-i o 

in 

r' 

rH 

o 




rH 

J. 


u 3 




n 

rj Eh 

Cl 

r-H 

t -J r H tN ro 




o 






t: 

C'r' 


• 

• 

o 

r j a; 

C 

o 

V£) 

O 

O 

ro 

o 


o 

O 

c 

t: 


vO 


‘G 

03 

V-t 

Q 03 

w 

H-» 


4 j 

4-3 


rH 


4 -> 

■iJ 

4 J 

t^- 

•H 

c 


si;- 

•J 

AJ 

tJ u 

u 


VO 









H 



U 

H 

rH 


l 4 

•rH 



ro 

CN 




U» 

r-j 

tn 


n rH o 'M "j 0? 

iJ d -H i> h 

wv > o u ^ 

fd -iJ ty ;vv,^ 
to ';} 0(4j 
c M c) ij: ^ 

n CO c 'HO 

O tj C'3 'D rd 


ti 

'■M 03 r3 

>.x: 

O 4J 4J T5 
< 'H O 
c/) ij r>i i-3 

c: 


O 4 -> 

r-i u 


o 


ra 


3 'O 
CT* 03 4v> 
O C IQ 
(U -r-t 03 
c > fd 
fd H 4J 

C >j 


u o 
03 o 
Oi c 03 

'i vU 

O' »d 
Es x: 133 d 
^ -H 'H 


C71 AJ [1 
t; 03 




( 3 ) 

0 T) 
E< C 

: j ti3 o 

o o 

I ic'C U3 
< U) 

1 < a: o 

< w 

H (; 

W H 

n 


v£) m CO v£) v£j <X) 
I CM 
V£> 


I — i r--( ) — i 

o o 


o ^ 

•U \£> 

I 

fO 


o 

u 


i O 

I 


rH G> lO 

r-f m 


f-H CD ro 

,-( I m 
( i 



■H 

•H Eh "O 

O H 

Q 

O 

rH 

^£^ VD 

O 

rH 






m 

ro 


rH 


£33 

fd is 


f-H is 

o 

rn 

ro 

r—i rH 

rT3 

m 










u 

O 

D *H 

03 

rH 

M D 











c 

D 


H 

o 

a 

< 

F, 

1h ^ 

£4 U 

c 

G 

H rH 

C 

u 






G 

03 


kX 

a 


£33 E < 

H 

<XJ t-H 

£0 .C 

n 

03 

0» G 

n 

03 






F) 

Q 


< 

03 

fd 

?d £X 

aJ 

c q: 

0( ,P 

h-3 

a 

< D 


Q 












Cx3 


Cl3 


















O 


U 

is 















r-«v 




z; 

O 











rH 

rH 

m 

rH 

u 


ID 


L3 

M 













rH 







F< a 

r-H 

rH (N r-' m CN 

rH 

rH 

O 









IG 


2: 


Eh 

o o 



(N • rH rH CN 



rH 





q 

r i 

c 

c 



o 


03 

U H 

q 

r-H 


c 

•H 






Id 

P 

p 

Id 

I 


?H 


CO 

Q u~t 

<d 

n 

H H c rH rH 

Id 

P 

H 





n 

F3 

F) 

ip 

• • 


Ft 



tiJ 



Id O. S3 G G 


D 

a, 









£33 Eh 


U 


23 

d: 



s: < F> F) H) 



< 







4J 


c < 


Zf 


u 











to 



kX 


O H 


a 


c 




iJ A-l 4-3 4J vU 



4J 




is 



03 


•H to 


U3 




< 


Qt Gi Gi Qi Gj 



>0i 




IH 



U 


4J ^ 


a 


w 


M 


O 03 C 03 03 



U 







X 


•H 





23 

rJ 


U C3 u u u 



u 




cu 



vD 


c: 





o 

< 


X X X X X 



X 




a 





'H 





rn 



03 03 03 03 03 

►H 


03 

E-h 

< 

ss 

X 




2 

UH 






Eh 



(S3 



CU 

M 

< 


u 



M 

03 





< 

CO 







Q 

a. 

4 

CO 



< 

Q 





Eh 

iD 






o 

;2 

< 

£h 

CiJ 



tx 






CO 

< 



rD 



CO 







CO 


unf.j 


PAGE 2b 


6. COMMUNICATIONS 

The communications center provides voice and teletype radio 
links to the SAO laser field stations in Arequipa and Natal. SAO 
also maintains FTS service within the continental United States, 
and is connected with Western Union, TELEX, AUTODIN, NASCOM and 
recently acquired RCA GLOBOM circuits for worldwide 
communications. The RCA GLOBOM circuit has been installed and 
operating on line since early December. We are presently sending 
75-80% of our overseas traffic, formerly passed by the W. U. 
TELEX, because of the less expensive message rates. W. U. TELEX 
is retained, however, for domestic use and some overseas traffic. 

The HF communication links with Natal and Arequipa have been 
operating satisfactorily; however, the teletype equipment 
continually requires adjustments and continues to be the weak 
link in the system. 

The SAO Cambridge radio room has undergone some major 
changes, including; the removal of the lOkw transmitter and 
associated power, the installation of a power distribution panel 
and circuit breakers for the remaining two Ikw transmitters and 
the receiver rack, the addition of an equipment interface 
switching panel and associated interface wiring, and the 
installation of a space HAL ST-6 demodulator/tone keyer. In 
addition to this, the equipment within the room has been 
rearranged to allow for more efficient operating procedures, and 
a new operations manual has been written and is in its final 
draft. A Bird model 4381 R.F. Power Analyst has been acquired 
for testing transmitter power at SAO Cambridge. 
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7. DATA SERVICES AND PROGRAMMING 


The Data Services and Programming groups perform central 
data processing and develop and maintain software for the laser 
tracking network. The Data Services group focusses on a 
continuous data exchange with laser sites in the SAO network and 
with other agencies and organizations both domestic and foreign. 
It examines raw ranging data, generates ephemerides for all 
satellites tracked by the network, and deposits the final SAO 
ranging data in the NSSDC. The Programming group is responsible 
for software used at all the laser sites in the SAO network, and 
at the central facility in Cambridge. 

7.1 Data serviges 

Data Services is responsible for two major areas of 
activity: the quick-look cycle and final data processing. On a 
weekly schedule, the SAO and cooperating laser stations send 
small subsets of their acquired quick-look data to Cambridge. 
These data form the basis for updated orbital elements, which are 
generated by Data Services and communicated back to the field 
stations for use in computing the predictions necessary for laser 
satellite ranging. 

The full data sets are recorded on Line magnetic tape and 
mailed from the SAO laser stations to Cambridge and sent through 
the final data-processing cycle. During this phase, data quality 
is assessed by an engineering filter, and the data are then 
processed through a noise filter, a time-corrections program, and 
a formatter, which outputs the data in the NASA coded format as 
well as the SAO coded format. It is the NASA coded formau that 
is supplied to the NSSDC. 
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The quick-look functions of the Data Services operation have 
evolved into a stable, reliable, and smoothly running procedure. 
Acquisition orbits were computed and transmitted each week 
virtually without incident. The quality of these orbits remains 
very high; ephemerides are now routinely computed to the 
sub-lC-meter level and, in the case of LAGEOS, to the 2-meter 
level. 


In the last six months, the Data Services group processed 
56,282 laser quick-look data points and handled 3,341 passes from 
the SAG and cooperating stations on GEOS 1 and 3, Starlette, 
BE-C, and LAGEOS. (See Table 1). 

A highlight of this reporting period was the participation 
of SAO in the Short MERIT Campaign of international 
investigations involving the Earth's rotation and pole position. 
The Data Services group, using Lageos data from the SAO and NASA 
laser networks as well as from certain cooperating foreign 
organizations, was able to obtain the earth's pole position as a 
by-product of the routine orbital determination starting in 
August, with the beginning of the 3-month MERIT Campaign. These 
pole positions were transmitted routinely to the B.I.H. in Paris 
as a rapid service to the world scientific community. 

Final data from the SAO laser network for the period March - 
April 1980 were transmitted to the NSSDC in July. May data were 
submitted in August, June data in September, and July 1980 
through October 1980 data were transmitted in December. 
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7.2 PI oa ramming support 

A small staff of computer programmers supports the operation 
of the laser tracking program. In addition to regular 
maintenance and upgrading of the minicomputer and production 
processing programs, the programmers develop software to meet new 
needs and support Data Services in routine processing as 
necessary. This group also analyzes test data for laser system 
maintenance and for planning laser system modifications to 
improve performance. The programming staff completed the data 
control and accounting package for the final data-processing 
cycle. This package facilitates archiving, control, retrieval 
and trouble-shooting final data. Basically this software 
maintains a file that allows the following of individual passes 
through the various processes, the observation of their status, 
and the generatation of statistics. This package was completed 
and placed in use by the Data Service Group. 

Our programmers completed work on the zero-set check of the 
orbit theory used in Cambridge (GRIPE) with that used at the 
field stations (FLPPS) . As noted in earlier reports, this is a 
full consistency check of quick-look and prediction software. At 
the present time an average of less than 10-arcsec directional 
difference and less than one hundred meters range difference 
between the full and abbreviated orbit packages has been 
attained. At the moment this level is adequate and the zero set 
procedure will be performed once a year or whenever changes are 
made to either package to ensure system integrity. 

However, in order to succeed in daylight ranging to LAGEOS 
it will be necessary to upgrade the orbit software and field 
predictions to the sub-microsecond level (100m) . In this 
semi-annual reporting period extensive investigation and 
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numericdi testing were performed on the compatibility between the 
headquarters orbital routines and the field programs. 
Discrepancies at this interface would be amplified with time over 
the prediction period. Problem areas that were identified 
included: 1), a small (10m) bug found in the Keplerian motion, 
2) , A 100-200m per week discrepancy found in the lunar theory 
used in the field prediction packages. Updating the field 
software will require a major programming effort, 3) , A 20m 
problem was found to exist in the tesseral harmonics. The most 
critical item is the lunar theory incompatibility, A major 
effort is now being applied to this issue. The other two 
problems which appear to be software "bugs" are being left in 
abayence at the moment. 

Having identified the major discrepancy as the lunar 
perturbation theory, it was decided to attempt to incorporate the 
fully developed integration theroy as used in the precision orbit 
determination program in the field minicomputer laser pointing 
software package (FLPPS), The technical problem will likely 
concern fitting this code into the limited memory available, and 
work is focussed in this direction at the moment. 

The other major area of intense programming effort concerned 
the laser control program which supplies the laser mount with 
prediction information on demand and logs the results of laser 
ranging on magnetic LINC tape for further processing. 

Two separate, but related hardware changes are in the works 
which requires substantial modification to the Direct Connect 
package. The first change is the replacement of the WD2000 
waveform digitizer with alternative waveform discrimination 
circuits. Considerable effort was made to develop a version of 
Direct Connect to properly treat this new data configuration. 


PAGE 30 


Preliminary versions for testing engineering were completed, but 
it is anticipated that continuing maintenance and modification 
will be needed. 

The second major change is the speed-up of the repetition 
rate to 30 pps. A stripped down version of Direct Connect was 
successfully tested in a simulator package (not an actual laser 
system) at the year's end. The next stages of development of 
this package will evolve by adding as many features as possible 
while still accommodating 30ppm. One casualty of this trade-off 
between processing speed and some initial "on-line" data 
processing was the second or reduced data file. The output of 
the Direct Connect program consists of a single LINC tape file, 
which will be a combination of raw and reduced data. 

The Quick-look Data Catalog program was completed this 
period, along with numerous minor modifications to field 
minicomputer and headquarters main-frame software. 
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8. DATA QUALITY 

Since the installation of the pulse chopper coupled with 
improved pulse processing, the SAO laser systems have been 
operating at an estimated 10 cm accuracy. This higher accuracy, 
however, must be monitored very carefully to ensure that the 
system does not degrade and that the ranging capability is 
maintained. In the following subsections, several key factors 
and test results associated with the determination of data 
quality are received Where possible, tests during this reporting 
period have been used to characterize the system at each station. 

8.1 gy.gtemfltig EuLOJ-a 

We routinely separate the systematic errors of the laser 
ranging system into three categories: spatial, temporal, and 
signal-strength variations. Spatial variations refer to 
differences in time of flight depending on the position of the 
target within the laser beam. Temporal variations relate to 
system drift between prepass calibration, satellite ranging, and 
postpass calibration. Variations in range due to changes in 
signal strength from pulse to pulse are a function of receiver 
characteristics and digitizer sampling interval. 

Spatial variations. Spatial variations, or the waveform 
error, were measured at Mt. Hopkins several times. Table 6 
gives the results before and after installation of the pulse 
chopper. The two wavefront measurements performed with the 
chopper show a maximum deviation of plus or minus 0.3 nsec (4.5 
cm) from the mean value across the wavefront. The standard 
deviation of the excursions is about 0.2 nsec (3 cm). 
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Table 6. 

Wavefront distortion. 




Spacing 

Average 

rms 


Date 


between 

number of 

wavefront 

Maximum 

of 


points 

photoelectrons 

distortion excursion 

measurement 


( arcmin) 

received 

( nsec) 

(nsec) 

26 Feb 1974 


0.3 

86 

1.4 

3.9 

18 March 1974 

0.6 

56 

O 

• 

00 

2.7 

9 NOV 197 8 


0.3 

56 

0.19 

0.6 

1 Dec 1978 


0.42 

87 

0.17 

0.6 

Temporal 

variations. Temporal variations are 

estimated by 

the difference between prepass 

and postpass 

calibration 

measurements 

' • 

Recent 

results from the 

three operating stations 

are shown 

in Figures 1, 2 and 3. 

These prepass 

and postpass 

calibration 

difference 

s represent an 

upper bound. 

since other 

statistical 

errors are 

also included. 
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Figure 1, 
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Signal-strength variations. Variations in apparent range 
with signal strength have beer examined with extended target 
calibrations over the dynamic raiige of the laser instrument. 
Figures 4, 5, and 6 illustrate recent calibration data taken in 
the high-satellite configuration. On high satellites (LAGEOS) , 
the system is operated right down to the single photoelectron 
level and the digitizer is operated with a 20-nsec baseline 
(1-nsec sampling) . In ) nis mode of operation, the extended 
target calibrations typically show a flat response in the 1 to 50 
photoelectron region to plus or minus 0.4 nsec (6 cm) or better. 
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Error budget. Usirg systematic error values of 4.5, 6, and 
4.5 cm for spatial, temporal, and signal-strength variations, 
respectively, for lew satellites and assuming that these errors 
are independent, the root-sum square (rss) error due to 
systematic sources is less than 9 cm (see Table 7) . For the 
higher satellites, the errors are typically 4.5, 6, and 6 cm, 
respectively, giving an rss of 9.5 cm. These, then, are the 
systematic errors that can be expected for data averaged over a 
pass. In addition, errors in timing, refraction, and spacecraft 
center-of-mass corrections need to be applied to the data; these 
errors are basically the same as those in the NASA systems, since 
essentially the same equipment types are used. 

Table 7. Summary of systematic errors. 



Low 

High 

Source 

satellites 

satell 


(cm) 

( cm) 

Wavefront distortion (spatial) 

4.5 

4.5 

System drift (temporal) 

6.0 

6.0 

Calibration (signal strength) 




rss 


8.7 


9.5 
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6.2 2y£i£ni Noise 

The range noise in the SAG lasers is determined in large 
part by return signal strength. Noise figures are typically in 
the 10- to 20-cm range for low orbiting satellites. With LAGEOS, 
at signal strengths of 1- to 5~photoelectrons, noise figures are 
several times larger (see Table 8) . It should be noted that the 
theoretical random statistics from photoquantization with the 
6-nsec chopped pulse are about 30 to 40 cm. This of course is 
degraded by jitter in the start and stop channel of the system 
electronics and by the current PMT. 
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Tabl 

LAGEOS passes from Arequipa, Natal 

Pass 

Stability 
( Pre-Post 
Calibration 


Date 

Station 

(CM) 

10/1 

Peru 

0.60 

10/1 

Peru 

1.50 

10/1 

Australia 

1.20 

10/1 

Australia 

-2.25 

10/2 

Peru 

-6.7 5 

10/3 

Australia 

-5.55 

10/3 

Peru 

6 .72 

10/4 

Australia 

-12.60 

10/5 

Australia 

0.60 

10/5 

Australia 

-2.25 

10/5 

Australia 

-3.15 

10/7 

Australia 

-4.05 

10/7 

Australia 

4.20 

10/7 

Peru 

-1.20 

10/7 

Brazil 

-7.05 

10/8 

Australia 

2.25 

10/8 

Peru 

5.55 

10/10 

Australia 

13.80 

10/10 

Australia 

2.10 

10/11 

Australia 

0.60 

10/11 

Peru 

-5.25 

10/11 

Brazil 

-2.25 


8 . 

and Orroral Valley, October 1980. 


RMS 


(Residual to 
Polynomial 
fit) 

(CM) 

Number 

of 

Points 

28.7 

37 

38.5 

52 

30.2 

63 

30.1 

169 

32.1 

90 

24.9 

34 

44.8 

51 

40.1 

47 

39.3 

13 

43.4 

98 

37.2 

113 

33.2 

128 

27.9 

34 

64.1 

7 

38.9 

13 

54.3 

47 

36.9 

45 

47.9 

75 

42.0 

66 

46.0 

36 

43.9 

42 

61.2 

15 
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Table 8 (Continued) 

LAGE03 passes from Arequipa, Natal, and Orroral Valley, October 1980. 


Pass RMS 

Stability (Residual to 
(Pre-Post Polynomial Number 


Date 

Station 

Calibration 

(CM) 

fit) 

(CM) 

of 

Poin 

lC/12 

Australia 

6.00 

26.5 

63 

10/12 

Australia 

5.55 

38.8 

74 

10/12 

Peru 

-4.35 

45.6 

31 

10/12 

Peru 

-5.25 

38.2 

92 

10/13 

Australia 

-4.05 

56.6 

62 

10/13 

Australia 

5.25 

27.0 

176 

10/13 

Peru 

3.60 

35.7 

11 

10/14 

Australia 

8.85 

57.0 

18 

10/14 

Australia 

6.00 

59.0 

52 

10/14 

Australia 

1.05 

39.4 

111 

10/14 

Peru 

-2.40 

59.0 

9 

10/14 

Brazil 

7.05 

50,0 

32 

10/15 

Australia 

-4.05 

36.0 

22 

10/15 

Australia 

-1,65 

45.40 

114 

10/15 

Peru 

-9.90 

36.5 

43 

10/16 

Brazil 

-10.65 

34.1 

9 

10/17 

Australia 

-8.10 

49.9 

82 

10/17 

Peru 

-7.20 

27.9 

63 

10/18 

Peru 

-0.60 

38.1 

88 

10/19 

Australia 

-0.60 

33.4 

63 

10/20 

Australia 

0.45 

36.1 

68 

10/20 

Australia 

-3.45 

28.2 

147 


Table 8 (Continued) 
LAGEOS passes from Arequipa, Natal, and Ocr 


Date 

Pass 

Stability 
(Pre-Post 
Calibration 
Station (CM) 

RMS 

(ResiduaJ 

Polynom; 

fit) 

(CM) 

10/20 

Peru 

1.05 

31.8 

10/22 

Australia 

3.45 

28.5 

10/22 

Australia 

-5.40 

31.5 

10/24 

Brazil 

-2.70 

45.6 

10/25 

Australia 

-0.30 

56.4 

10/25 

Australia 

9.60 

50.0 

10/26 

Australia 

-5.55 

23.0 

10/26 

Australia 

-0.7 5 

39.7 

10/26 

Australia 

4.05 

65.5 

10/28 

Australia 

7.05 

52.9 

10/28 

Australia 

7.35 

38.9 

10/29 

Australia 

-7.65 

28.4 

10/29 

Australia 

6.75 

36.5 

10/29 

Australia 

-3.00 

33.8 

10/29 

Brazil 

-4.95 

37.3 

10/30 

Australia 

7.05 

30.2 

10/30 

Peru 

-2.40 

30.1 

10/31 

Peru 

0.00 

51.5 
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8.3 Lffiaa Asm Mialysla 

In addition to the system calibration tests and short arc 
analysis, we also examine the data with a long arc analysis to 
determine an upper oound on systematic bias. Using a 10-12 day 
arc, we determine a mean offset for each pass which is the result 
of instrumentation errors and modelling error due primarily to 
gravity field and station position. Two such runs during this 
reporting period are shown in Tables 9 and 10 Most of the pass 
offsets denoted as <p> are in the neighborhood of 10-20 cm or 
less, with some examples as large as 50 cm. At the moment this 
type of behavior is anticipated from the modelling errors above, 
so the current feeling at both GSFC and SAG is that the ground 
based measurements are indicative of the ranging accuracy. 
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Table 9 
LAGEOS 

Long Arc Analysis 
Final Data 

September 12-25, 1979 


STATION 

DATE 

HOUR 

NO. PTS. 

SIGMA 

<P> 

7 907 

Sept. 13 

01 

30 

48 cm 

9 cm 

7907 

Sept. 13 

23 

46 

47 cm 

- 2 cm 

7907 

Sept. 14 

09 

21 

4 8 cm 

-53 cm 

7907 

Sept. 16 

10 

15 

66 cm 

-41 cm 

7907 

Sept. 19 

00 

34 

47 cm 

"14 cm 

7 907 

Sept, 19 

09 

24 

42 cm 

-13 cm 

7907 

Sept. 21 

23 

42 

47 cm 

5 cm 

7907 

Sept. 22 

08 

19 

67 cm 

- 1 cm 

7907 

Sept. 24 

00 

24 

51 cm 

-11 cm 

7907 

Sept. 24 

23 

8 

48 cm 

4 cm 

7 943 

Sept. 12 

10 

48 

31 cm 

- 3 cm 

7943 

Sept. 15 

19 

38 

55 cm 

-11 cm 

7943 

Sept. 18 

09 

51 

24 cm 

-17 cm 

7943 

Sept. 23 

09 

35 

32 cm 

- 6 cm 

7943 

Sept. 25 

19 

26 

37 cm 

30 cm 


7907 Arequipa 
7943 Orroral Valley 
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Table 10 


LAGEOS 

Long Arc Analysis 
Final Data 
August 4-14, 1979 


STATION 

DATE 

HOUR 

NO. PTS. 

SIGMA 

<p> 

7091 

Aug . 

7 

02 

29 

51 

cm 

-17 

cm 

7082 

Aug . 

7 

05 

31 

27 

cm 

21 

cm 

7082 

Aug . 

10 

04 

38 

6 

cm 

4 

cm 

7082 

Aug . 

11 

08 

30 

30 

cm 

32 

cm 

7082 

Aug . 

14 

03 

31 

14 

cm 

- 9 

cm 

7943 

Aug . 

4 

10 

42 

22 

Cm 

- 4 

cm 

7943 

Aug . 

7 

10 

42 

30 

cm 

3 

cm 

7943 

Aug . 

13 

08 

44 

37 

cm 

6 

cm 

7943 

Aug . 

14 

20 

21 

40 

cm 

5 

cm 

7907 

Aug . 

4 

03 

13 

43 

cm 

1 

cm 

7 907 

Aug . 

5 

01 

21 

29 

cm 

- 5 

cm 

7907 

Aug . 

8 

00 

13 

45 

cm 

- 3 

cm 

7907 

Aug . 

8 

23 

12 

45 

cm 

-34 

cm 

7907 

Aug . 

11 

00 

4 

18 

cm 

-21 

cm 

7907 

Aug . 

11 

23 

5 

39 

cm 

•• 2 

cm 

7 907 

Aug . 

13 

01 

8 

39 

cm 

9 

cm 

7907 

Aug . 

13 

23 

8 

33 

cm 

- 3 

cm 

7 92 9 

Aug . 

11 

22 

12 

34 

cm 

-29 

cm 


7091 Haystack 
7082 Bear Lake 
7943 Orroral Valley 
7907 Arequipa 
7929 Natal 
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9. RELOCATION OF THE MT. HOPKINS LASER 

The NASA Crustal Dynamics project calls for the development 
of a system of mobile lasers (MOBLAS) to be used in conjunction 
with lasers already in place to achieve a global network. Early 
plans included mobile lasers in western Australia, Diego Ga’’cia, 
Kwajalein, American Samoa, Marcus Island, and Madrid, Spain; the 
SAO laser sites at Arequipa, Natal, and Orroral Valley; and the 
University of Hawaii site at Mt. Haleakala. This development 
plan has since been modified to use the Institut fur Angewandte 
Geodasie's laser in Wettzell for coverage over Europe, thereby 
obviating the need for a mobile laser in Madrid, to delete Marcus 
Island, and to relocate SAO's Mt. Hopkins laser to India, rather 
than placing a mobile laser at Diego Garcia. This last decision 
was based on the poor climatic conditions at Diego Garcia, the 
inability to monitor local deformations there, and the high costs 
of operating a system on the island. 

In our FY 1980 proposal, SAO outlined a program to relocate 
the Mt. Hopkins laser station to a site in India. The plan was 
based on a cooperative arrangement with the Indian Space Research 
Organization (ISRO) at a mutually agreeable site. SAO proposed 
to provide Excess Currency (PL480) to partially support the 
operation. ISRO suggested Kavalur, a site near Bangalore where a 
Soviet laser and camera are presently being maintained. After 
consideration of logistics, facilities, weather, and location it 
became evident that Kavalur would not be an adequatre site and 
SAO proposed that the laser be located at the Uttar Pradesh State 
Observatory at Naini Tal, a site which has been an SAO Baker-Nunn 
camera site for over 20 years. The site was deemed acceptable to 
SAO because of its excellent weather conditions, its established 
facilities and personnel, and its accessibility in terms of 
logistics and communications. Earlier concerns about the 
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degrading effects of local crustal motions on long baseline 
measurements have been allayed with a plan to use geodetic and 
geophysical measurements to monitor local crustal movements. A 
proposal was submitted to the Smithsonian Institution for PL-480 
Excess Currency to support part of the expenses to be incurred by 
SAG, and plans were developed for the required modification of 
the Baker Nunn building at Naini Tal for laser installation. 

ISRO did not provide information on any candidate laser 
station site other than Kavalur; in mid-February 1980, NASA and 
SAG agreed that the UPSG site at Naini Tal would be suitable for 
the laser station and communicated this decision to ISRG and 
UPSG. A revised MGU based on placing the laser at Naini Tal was 
drafted by SAG and sent to ISRO for comments. The revised MOU 
was discussed at a meeting with ISRO at Bangalore, India, in May 
1980. It appeared that there was general agreement on the roles 
that ISRO, UPSG, and SAG would play in the cooperative laser 
program. It was anticipated that all participants would give 
final approval by late summer. 

Under the arrangements discussed with ISRO, UPSG and the 
U.S. Embassy, it was anticipated that: 

1 . ISRO would pay the bulk of the local expenses to be 
incurred at Naini Tal. 

2. SAG through Excess Currency Program (PL 480) would 
support transportation of goods, travel expenses, local 
administration, and some communication and other local expenses 
in India. 

3. SAG through NASA funding would pay salaries and some 
expenses for U.S. employees, plus station support costs normally 
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incurred within the United States. 

4 . The U.S. Embassy in New Delhi, wculd provide the 
necessary administrative and communications support to bring this 
project to fruition. 

In the interim period, there has been some very limited 
communication with ISRO, mainly through the Embassy, which 
continues to convey a favorable position for ISRO. Discussions 
between NASA and ISRO, including discussion between Drs. Frosch 
and Dhawan, indicate that ISRO is still giving this program very 
serious consideration and a response will be forth coming in 
early 1981, after they have completed their long range planning. 

Assuming an approval for the program in January, a revised 
schedule for the establishment of the laser station in Naini Tal 
is shown in Figure 7. The schedule is given in terms of calendar 
months. Based on these assumptions, our best estimate at this 
time is that the India station will be in operation by early FY 
1982. 
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10. LASER UPGRADING 


Under revision 1 to the Fyi981 Proposal (P987-6-80) and 
subsequent approval by NASA, SAG has embarked on a program ol 
laser upgrading to bring the SAG lasers co 3-5cm accuracy and to 
improve data yield for the crustal dynamics program. The key 
elements of the program are to reduce the laser pulse width from 
6 nsec to 2-3 nsec by modifying the pulse chopper and increasing 
the pulse repetition from 8ppm to 30ppm to increase the data 
yield. To accommodate the increased pulse repetition rate and to 
reduce range noise, we are replacing the waveform digitizer with 
an analog pulse processor and by using a better PMT and PMT base. 
System Signal to Noise, particularly during twilight and daylight 
ranging, is being increased with a redesign of the photoreceiver 
to accommodate a 3 Angstrom interference filter and a 
modification to the range gate generator to reduce the gate 
window by a factor of 10. Associated into these changes are 
modifications to hardware, software, and operational procedures. 

The modifications are being implemented primarily for 
LAGEOS. All of the modifications will be applicable to some 
extent to the lower satellites; however, the higher angular 
tracking rates and the greater uncertainty in orbital position of 
the lower satellites will preclude some of the potential for 
higher data rates and narrow range gate windows. 

The program will take approximately 18 months, and should 
have the lasers prepar-d for the next MERIT Campaign. 
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10.1 Increased Data Rate 


The major limitations in pulse repetition rate for the SAG 
lasers has been the mount and the pulse processing in the 
Computer. In the case of LAGEOS, where the angular rates are 
very slew, the mount can operated at rates well in excess of 30 
ppm and the computer becomes the limitation. The implementation 
of the analog pulse processor ameliorates this problem. 

Increasing the laser pulse repetition rate from 8 ppm to 30 
ppm requires modifications to (1) the laser power supply and (2) 
the laser control unit. All other hardware is presently 
compatible with the increased rate. The modifications to the 
power supply are being done by the SAO engineering staff with 
modules and components to be furnished by Laser Industries, Inc. 
Technical discussions are presently underway with Laser 
Industries. The current plan is to have them under contract in 
early CY 1981 for the required hardware. The modifications to 
the laser control units require that these units be replaced with 
new laser control units. These are now being designed and will 
be built under the direction of the Senior Electronics Engineer. 


A prototype software package which runs at 30 ppm with 
simulated analog pulse processor data has been written and is 
operating at SAO. Efforts are underway to add operational 
features to the program. Changes in format and some control 
considerations will make it necessary to change some of the 
subsequent data processing software once the field program are 
completed . 
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10.2 Reduction in Laser Pulse Width 


The 
building 
chopper . 
builders 
Blumlein 


laser pulse width will be reduced to 2-3 nsec by 
a shorter ceramic Blumlein structure for the pulse 
Discussions are underway with Lasermetr ics , Inc., the 
of SAO pulse choppers, who will fabricate the new 
structure . 


10.3 Improved Range Gate 

The SAO lasers currently use a range window of 6-10 
microseconds. This poses no difficulty during the nighttime 
hours; however, at twilight and under daylight conditions the 
system is very noisy, and operation on LAGEOS is precluded. To 
improve the system, SAO is modifying the range gate generator to 
oprated at windows down to 0.1 microseconds. Design of this 
modification to now underway. Some reformatting in the hardware 
and the addition of som circuit components. 

In addition to the hardware modifications, the prediction 
software is also being modified to be compatible with 
submicrosecond predictions (see Section 5) . 


10.4 Modifications to the Photoreceiver 

Improvements are also being made to the photoreceiver to 
improve range accuracy and the return signal to noise ratio. A 
new PMT (Amperex 2233) and a new PMT base have given considerable 
improvement in pulse shape and repreducibility . In addition, the 
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photoreceiver is being modified to accommodate a fast pulse 
shutter and a 3 Angstrom narrow band filter. This should provide 
an additional factor of 2.6 improvement in signal to noise from 
the present 8 Angstrom filter, and will offer protection for the 
PMT from daylight and laser backscatter ing . 

The design of photoreceiver modifications has been 
completed. Vendors have been selected for the components to be 
ordered in January. The 3 Angstrom filter will be purchased from 
Daystar Filter Corp. The shutters, already tested at Mt. 
Hopkins, will be purchased from Vincent Associates. 


10,5 Analog Pulse Detector 

A fundamental element to the increased data rate and the 
improved accuracy is the analog pulse detector. Not only does 
this reduce considerable the data processing per pulse over the 
digitizer and thereby make the 30ppm mode possible, but for very 
narrow pulses it also avoids errors introduced through inadequate 
pulse amplitude sampling. In addition, it also avoids the 
inherent complexities of the WD2000 and the related pulse 
processing software. 

SAG has built three different analog pulse processors which 
have been undergoing testing at Mt. Hopkins. All appear to give 
results at least recomparable to the WD2000; detailed tests are 
underway at Mt. Hopkins to demonstrate if one of the systems, or 
combination thereof is superior. The first processor uses a 
match filter (matched to the laser output pulse) with 
differentiator and a cross-over detector (discriminator with a 
low threshhold and a slope dependent arming circuit) . The second 


PAGE 56 


processor is a tapped delay line integrator with a cross-over 
detector. With this scheme, the taps from the delay line are 
summed and differenced to give effectively a pulse area 
differentiation which is then examined for the point of 
cross-over. The third detector is based on constant fraction 
discriminators. 

Extended target calibrations are being used to assess the 
performance of each system, and to make a final decision on the 
system configuration for both the start and stop channels. We 
anticipate that the tests will be completed by late January and 
that a decision on final system configuration can be made in 
February after all of the data are evaluated. Many of the parts 
for the production systems have already been ordered. 
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11. PERSONNEL 

11.1 Visitors 

Several companies, who are potential Service Contractors in 
support of the SAO Mission Contract bid, were represented here at 
headquarters as SAO commenced meetings for the discussion of the 
Mission Contract proposal. 

In August, Memphis Norman, a budget examiner from the Office 
of Management and Budgeting, came to Cambridge for general 
administrative discussions. 

A number of field personnel visited Cambridge during this 
reporting period. Stephen Criswell an Donald Patterson visited 
during August while on home leave. Dana Seaman arrived in 
September for several days of technical discussions. William 
Dennison Rich, manager of the Brazil station, and David 
Hallenbeck, manager of the Peru station, took annual leave during 
November, and also visited SAO headquarters for discussions. 


11.2 F £LS Q nnfi.l 

In August, Mr. Dennis Anttila was hired as a new Observer 
for the Arequipa station; however he resigned shortly 
thereafter. 

There were a few changes in headquarters personnel during 
this reporting period. In August, Eileen Brady joined our staff 
as a secretary. David Kirk Gilmore transferred from his position 
of Computer Technician in STADAD to the Optical and infrared 
Astronomy Division at Mt. Hopkins, where he is a Physical 
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Science Technician. Jean Kintner joined our staff during 
November as a computer programmer, transfering from another part 
of the Observatory. Suzanne Romaine also transfered to another 
division, leaving the position of a computer programmer; Saeqa 
Oil replaced Suzanne during the latter part of December. 

Walter Munn, Master Observer at Wallops Island, Virginia*, 
retired at the beginning of December. 


11.3 Tia v el 

Various headquarter personnel went to Perkin Elmer, an 
engineering firm in Wilton Connecticut, to deliver and pickup the 
photoreceiver mirrors undergoing modifications there. 

Dr. Michael Pearlman attended various conferences and 
participated in administrative and techrsical discussions at NASA 
headquarters and Goddard Space Flight Center. Margaret Warner 
visited Wallops Island and the Smithsonian Institution during 
this reporting period. 

Two of our engineers traveled to the Mt . Hopkins station 
during this six month period. James Maddox worked on upgrading 
of the photoreceiver and carried out tests with the PIN diodes 
during his stay in November. Noel Lanham arrived shortly 
thereafter for field trials with the prototype analog detector 
system and to perform maintenance on the Mt . Hopkins laser 
system in readiness for relocation to India. 


* 


Supported by NASA Contract NAS6-2910 
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Semiannual Report No. I 


1. INTRODUCTION 

This Semiannual Report covers the initial part of work done 
under NASA Grant NAG5-36, Study of the Feasibility of 
Satellite-to-Sateliite Doppler Tracking for Mapping of the 
Earth's Gravity Field Suitable for Geophysics and Oceanography, 
and covers the period 1 April 1980 through 30 September 1980. 
The overall objectives of this grant ate summarized in the next 
section. During this first reporting period, work was undertaken 
in support of three tasks. First, a general analytical study was 
undertaken to obtain simple relations about the sensitivity and 
resolution of the gravity field determination to three different 
types of observations. This study allows relatively inexpensive 
comparisons of different types of data and provides guidance for 
more detailed numerical modeling experiments. Second, careful 
conversion and upgrading of the computer software to generate 
simulated satellite to satellite tracking data. Third, initial 
development of the analytical solution for the satellite 
perturbations due to the anomalous geopotential. This analytical 
solution forms the basis of the Generalized Inversion Method to 
be used for mapping the geopotential at the earths surface from 
doppler measurements at satellite altitudes. 

In addition during this period, SAG presented preliminary 
results at a GRAVSAT working group meeting at Goddard Space 
Flight Center on March 12-13, 1980. 
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2. OBJECTIVES 

The Smithsonian Astrophysical Observatory (SAO) is 
supporting the NASA Geodynamics Program. One of the goals of 
this program is to use space techniques to study the fine 
structure of the earth's gravity field in length scales between 
100 and 2000 km. The gravity field (or external potential) 
provides one of the few observations related to the earth's 
internal mass distribution and, therefore, the internal structure 
of the earth. Equally important is that the geoid (a surface 
defined by the external potential) the the reference to which 
changes in sea-surface elevation are referred. These sea-surface 
elevation changes can be measured with sa tel lite-to-sea-surf ace 
altimeters such as Geos 3, Seasat, or Topex, and they contain 
information of the ocean circulation and eddy field. SAO will 
pursue the three topics summarized below: 

1. We will investigate analytical solutions for the 
relative velocity of such configurations as high-low and low-low 
satellite systems. The relative velocity of two satellites is 
determined by the anomalous potential of the earth, among other 
things. Given a description of the anomalous potential, we can 
calculate the relative velocity to be measure with a SSDT system 
(the forward problem). These analytical expressions will be used 
to study the sensitivity of the SSDT system to the anomalous 
potential as a funcf'on of amplitude, wavelength, and system 
conf iguration. 

2. Observations of relative velocity at satellite altitudes 
can be used to calcuate the anomalous potential at the earth's 
surface (the inverse problem). Three steps are required: 1) 
solution of the boundary value problem, 2) solution of the 
downward continuation problem, and 3) regularization at the 
earth's surface. Using analytical solutions for satellite 
motion, we propose to investigate two methods to solve these 
problems - generalized inverse methods and least squares 
collocations methods. 
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3. We propose to investigate a novel mission concept that 
uses two low satellites with moderate eccentricity. The two 
satellites would be in phase^ in the same piane^ but with their 
arguments of perigee differing by 180 degrees. Therefore, when 
one satellite is at perigee, the other will be at apogee, and the 
SSDT will measure the relative radial velocity (the high-low 
configuration). Approximately one-quarter revolution later, the 
SSDT will measure the relative along track velocity (the low-low 
configuration). Because of the secular motion of perigee, this 
system would sample the whole earth with measurements of both 
radial and along-track relative velocities. Such a mission 
conf igi'ation would also have a considerable advantage in terms of 
system lifetime and resolvability of the geopotential. We intend 


to carry out a 
including a study 
gcopotential . 


system simulation of 
of the recovery and 


this mission 
separability 


concept , 
of the 
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3. ANALYSIS 

The lack of simple closed analytic expressions relating the 
observable, dr/dt, with the anomalous potential, U, has hampered 
a careful discussion of a GRAVSAT mission. Consequently two 
schools of analysis have developed: a) those who plan and carry 
out large computer simulations (Douglas et.al. 1980, Schwarz, 
1970, 1972), and b) those who extend very simplified models that 
are tractable (Breakwell, 1979, Fischell, 1978). Results from 
these simplified models can give useful insites and provide some 
guidance in planning extensive simulations. They are also useful 
in developing mission strategies and profiles. Ultimately one 
expects these two approaches to agree, and then be able to select 
the particular technique appropriate for the problem at hand. 

The analysis described here follows the ideas of Breakwell 
(1979). The extension to variable satellite spacing, satellite 
configuration, non - polar orbits, mixing of data types, and 
averaging methods are given, with some examples. 

The satellite is assumed low, h/R << 1. , and the analysis 
is restricted to short wavelengths. In this case a flat earth 
approximation can be used and the estimation problem (downward 
continuation) is handled by two dimensional fourier transforms. 

One can begin with the upward continuation of the surface 
potential using the flat earth form of Poisson's integral 
formula: ^ ^ 

U(x,y,z2) * - I (z2 - zl)U (x,y,zl) * [x^ + y^ + (z^ - zD^J dA 


where the height is in the z direction. If the two dimensional 
fourier transform of the potential is U(wx,wy,z) the fourier 
transform cf (1) is found to be: 

U(wx,wy,z2) = exp(-(z2-zl)w) * U(wx,wy,zl) 

where w**2 =• wx**2 + wy**2. The tranforms for various 

observables are easily found (following Breakwell) , with 
U0=U (wx, wy , 0) and h=z2 as: 

Gradiometer 

(Uxx,Uyy,Uzz,Uxy,Uxz,Uyz) » 

U0*exp(-wh) * (-wx**2,-wy**2,-w**2 ,-wx*wy, -i*wx*w, -i*wy*w) 

Velocity measurements 

(Vx , Vy, Vz) =U0*exp(-w*h) * (1/V, wy/wx/v, -w/wx/V/i ) 

where i**2=-l. The possible measurements are combinations of 
these observables. The measurements considered are given in 
table 1. In these cases, the LO-LO is considered as two 
satellites at the same altitude (h) separated by a distance (A) 
traveling in the north south direction, i.e. parallel to the x 
axis. The HI-LO case is considered as two satellites, one above 
the other, at altitudes hi and h2. For the preliminary analysis 
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considered here, the eccentric LO-LO is considered as 
contributing part of the data in the pure LO-LO configuration, 
and part in the pure HI-LO configuration. Note, that the non 
polar configurations are obtained by taking the appropriate 
combination of the x and y components. 



PAGE 7 


The potential U(x,y,z) can be ised to calculate surface 
data, gravity anomalies or geoid heights. The fourier transform 
of these surface data can also be calculated by straightforward 
formulae. However, the quantities of interest are averages of 
gravity or geoid. The traditional averaging algorithm is to 
average a square block of side a. These blocks are easy to 
calculate but difficult to deal with theoretically. An 
alternative is to chose a spherical cap of radius a. This gives 
relatively simple theoretical formulae, but one rarely actually 
deals with data averages in this way. A third way, v;as chosen by 
Breakwell for its "simplicity". It is a Gaussian smoothing 
function. In fact all these methods are easily treated, within 
the framework outlined above. We give in table 2, the transfer 
functions for these averaging methods. 



Gravity Anomaly Transfer function 
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We are now in a position to estimate some scalar function 
from the random process UO by: 


'f (o) = f (u )U0 (u ,h) 
from a set of measurements by: 
Jf (u)) = El(u-)U0(o)+w 


( 2 ) 


(3) 


where H(u) is given in table 1, and f(u)> is given in table 2 , and 
W is the fourier transform of the measurement errors. If go 
the spectral density of the surface potciitial variations, and 
4i is a matrix of measurement error spectral densities, then an 
estimate of the spectral density of the estimated quantity 
(gravity anomalies or geoid heights) ^ is given by: 



-^UO (co) f f ( 


( 4 ) 


For this analysis we use Kaula's celebrated rule of thumb, 
UO (1 ) =GM/a *3 0* * -5/1 * *2 , for the spectral density of the surface 
potential . 


Equations 1, 2, and 3 can be used to compare different data 
types, and distribution. Given the assumptions, these formulae 
are unlikely to give absolute accuracies. However, they should 
provide some assessment of the relati’^e accuracy of some 
measurement systems, and some different spectral characteristics 
of these systems. 


The first comparison is of the various averaging methods, 
and the sensitivity of the measurement systems. These results 
assume a polar satellite, a data rate of one per 10 seconds, and 
a mission life time of 90 days. Since all the statistical 
parameters scale with these figures, one can modulate the 
absolute magnitude of the result very easily. However, the 
relative size of the errors is insensitive to changing these 
parameters. The measurement accuracies are given for 
conservative and optimistic estimates for both gradiometers and 
relative velocity measurements. For purposes of establishing 
nomenclature recall that one Eotvos Unit is 10**-9 
(cm/sec* *2 ) /cm, and that one micron/sec is 10*'*-3 cm/sec. Table 
3 gives the uncertainty of the gravity anomalies or geoid heights 
for a variety of averaging methods, and data types. Two points 
should be noted. First, the variety of averaging algorithms, 
give comparable results, and any one can be used for purposes of 
comparisons. Second, according to this simulation any of the 
proposed systems could achieve the mission objectives of a 
GRAVSAT. 
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The gradiometer and satellite to satellite tracking 
measurments sample the geopotential spectrum in different ways. 
Figures 1, 2, and J display this sensitivity for a number of 
different data accuricies and orbital configurations. Figure 1 
contrasts gradiometer accuracy, altitude, and sampling rate. The 
percentage of each harmonic of degree 1 that is recovered is 
plotted against degree. Clearly, increased accuracy, higher data 
rate, and lower altitude give greater accuracy. Modulating any 
of these mission parameters can improve (degrade) the accuracy. 
However, the shape of the curve is unchanged. Modulating these 
mission parameters only shifts the curve. 

In figure 2, the LO-LO configuration is examined for 
different accuracy, altitude, and separation. Again, varying the 
mission parameters can improve (degrade) the result. However, 
varying separation distorts the pectral response. In curves 3 
and 4, for 300 km spacing, the sensitivity is degraded. This is 
basically related to the well known fact that the two satellite 
system is insensitive to geopotential wavelengths comparable to 
the spacing of the satellites. 

Figure 3 shows the spectral response of the eccentric LO-LO 
system. To some extent the reduced sensitivity to wavelengths 
comparable with the satellite spaceing is somewhat mitigated by 
having the verticle velocity measurement. However, there is 
still some distortion of the frequency response due to this 
commensur ability . 

The frequency response of the three measurement systems can 
be contrasted in another way. Figures 4, 5, and 6 plot the error 
or uncertainty in each wavenumber (1) and a function of 
wavenumber, for the gradiometer, the LO-LO and the eccentric 
LO-LO. Here the frequency response is graphically displayed. 
Figure 5, shows of the uncertainty growing without limit for 
wavelengths corresponding to the satellite spacing. One would 
conclude that if the LO-LO system is adopted, then the spacing 
should be modulated to sample all wavenumbers equally. The 
gradiometer system has a somewhat more monotonic response. 
Especially for the higher degree gravity anomalies the 
gradiometer seems to have an advantage. The eccentric LO-LO 
shows some of the insensitivity at middle wavelengths as the 
baseline LO-LO system. However, the frequency response is 
broader, and finetuning the mission profile may result in a 
uniform spectrum. 

In figure 7 an overall comparison of the system response of 
all three measurement systems is given. With the adopted mission 
parameters, the gradiometer seems to have the best performance. 
However, all of these curves are relative. This figure is more 
useful to estimate the tradeoffs for spectral response, than for 
absolute accuracy estimates. For example, it is clear that when 
the gradiometer sensitivity begins to decrease, it fall more 
sharply than, say, the eccentric LO-LO system. 
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4 . SOFTWARE DEVELOPMENT 


Software development for the analysis can be divided into 
two categories. The first consists of the modification and 
testing of already existing routines. Under this heading fall 
the various routines that make up a numerical integration 
software package for calculation of satellite trajectories. 
Briefly, this program takes as input the initial Kepler elements 
or state vectors of one or more satellites and integrates the 
equations of motion to obtain the positions and velocities of the 
satellites at any later time. Output includes the relative 
line-of -s ight velocity between the satellites in the relevant 
case of two satellites. The force acting on each satellite 
includes the force of the earth's gravitational field (central 
and non-central parts) and various other effects that may be 
switched on or off as desired. These additional effects include 
radiation pressure, atmospheric drag, and the gravitational 
attraction of the sun and moon. An additional feature is the 
option of including surface layer density gravity anomalies of 
the earth's gravitational field of any dimension and intensity at 
any position on the earth's surface. This feature is of course 
particularly relevant to the present analysis since the detection 
of such gravitational anomalies is the object of the low-low 
Doppler experiment. The integrator program can be used to 
generate reference data, against which we can compare the results 
of the analytically determined values of the relative velocities, 
thus obtaining a measure of the accuracy of our analytical 
formulae. The integrator could also be used in conjunction with 
the inverse problem analysis. 

The integrator had originally been written for a CDC 6400 
computer. Most of the routines had been converted for use on the 
SAG VAX 11/780 computer, but some had not. The first requirement 
was to make sure all the routines were converted and that the 
program was working correctly. This has been accomplished. The 
critical quantity for our analysis is the relative line-of-sight 
velocity between the two satellites. To test whether the 
integrator was performing its job with sufficient accuracy, we 
used the accuracy of the relative velocity as the criterion. We 
compared relative velocity output from the integrator program 
with the analytical values obtained by solving the Kepler 
equation for corresponding times, given the same initial 
conditions, for the case of simple two-body motion. The accuracy 
was well within the limits required. Typically, the difference 
in relative velocity obtained by the two methods was 0.03-0.04 
micrometer/sec, corresponding to a relative error of 10'° to 10 '7. 


PAGE 21 


The next stage in testing the integrator routine v'as to 
determine the sensitivity of the integrator to graininess of the 
surface layer gravitational anomaly. For example, the effect of 
a single 110km x 110km surface anomaly on the the satellite orbit 
is not the same as that of the same surface area with the same 
intensity divided into 9 equal surface patches, and so on. What 
we needed to determine was the limit at which further division of 
the surface area into a smaller grid no longer caused a change in 
the relative line-of-sight velocity (the experimentally 
determined quantity in a Doppler experiment) greater than the 
hoped-for accuracy of the experiment. This meant the point at 
which the change brought about by a finer grid was less than a 
micrometer /sec . 

The results of this study were that with both satellites at 
200 km altitude a 100 mgal anomaly on a 110km x 110km surface 
patch required an 11x11 grid (121 subpatches) to obtain the 
desired accuracy in relative velocity, and at J50 km altitude a 
13x13 grid was required. Since the computer time required rises 
with the number of surface patches, this result emphasizes the 
desirability of obtaining analytical formulas for the relative 
velocity components, which is one of the primary goals of the 
study. 

The second category of software development is that of new 
computer programs to be used in the determination of relative 
velocity components as a function of time (and position) for an 
arbitrary surface anomaly distribution and for the inverse 
problem of determining the surface anomaly distribution from the 
relative velocity. At present, we are in the process of 
developing a program to determine the relative velocity 
components for the case of area mean anomalies using formulas we 
have derived. The first tests will be made with the simple 
two-body central force as the background force, with the 
important J2 harmonic term added later. configurations. 
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5. PLANS 


The program we plan to follow in order to achieve the objectives 
outlined above may be briefly summarized: 

1) Continue to develop software for obtaining perturbed relative 
velocity components with minimum expenditure of computer time in 
order to make computer simulations practical for general surface 
layer density gravity anomaly distributions as outlined above, 

2) Concurrently to develop and test inversion methods for 
obtaining gravity anomaly distributions from simulated velocity 
data, studying the sensitivity of the results to variations of 
relevant parameters such as satellite height, orbital parameters, 
and the intensity and distribution of density gravity anomalies. 

3) To precede by studying the simplest cases first, i.e., simple 
two-body motion plus mean surface layer anomalies, and then 
adding more perturbative terms to obtain progressively more 
realistic models as each level or approximation is understood for 
the different orbital configurations. 


I'AGK 2 3 


6. REFERENCES 


Bceakv/ell, J.V., Satellite 
Gravity Variations, 
Vol.xxvii, no. 4, pp. 


Determination of Short 
Journ. Astronautical 
329-344, Oct. -Dec. 1979 


Wavelength 
Sciences , 


Douglas, B.C., Goad, C.C., and Morrison, F.F., Determination of 
the Geopotential from Satellite-to-Satellite Tracking Data, 
Journ. Geophys. Res., vol.85, no. BIO, pp. 5471- >480 , 10 

Oct. 1980 


Fischell, R.E. and Pi.acane, V.L., A Drag-Free 
System for Improved Gravity Field 
Appl ications of Geodesy to Geodynamics, 
Mueller, Reports of the Department of 
Ohio St. University, No. 280, pp. 213-219, 


Lo-Lo Satellite 
Measurements in 
edited by I. I. 
Geodetic Science, 
1978 


Schwarz, C.R., Gravity Field Refinement by Satellite-to-Satellite 
Doppler Tracking, Dept. Geodetic Sci., Ohio State 


University , 

Res. Found. 

Rep. 

No. 

147, December 1970 

Schwarz , C. R. , 

Ref inement 

of 

the 

Gravity Field by 


Satellite-to-Satellite Doppler Tracking in The Use of 
Artificial Satellites for Geodesy, Geophysical Monograph 
15, American Geophysical Union, Washington, D.C., 
HenriKsen, S.W., Mancini, A., and Chovitz, B.H. editors, 
1972. 


APPLNDIX B 


INVESTIGATE THE TIME EVOLUTION OF THE CORRELATION BETWEEN BATHYMETRY 
AND THE GEOID HEIGHT IN THE PACIFIC OCEAN 

tjntract NA79SAC00739 
Quarterly Reports Nos. 4 & 5 
For the period 1 May through 31 October 1980 


Principal Investigator 
Dr. Micheline C. Roufosse 


Prepared for 

National Oceanic and Atmospheric Administration 
Rockville, Maryland 20852 


December 1980 


Smithsonian Institution 
Astrophysical Observatory 
Cambridge, Massachusetts 02138 


The Smithsonian Astrophysical Observatory 
and the Harvard College Observatory 
are members of the 
Center for Astrophysics 


INVESTIGATE THE TIME EVOLUTION OF THE CORRELATION BETWEEN BATHYMETRY 


AND THE GEOID HEIGHT IN THE PACIFIC OCEAN 
Quarterly Report Nos. 4 & 5 
Contract NA79SAC00729 

For the period 1 May through 31 October 1980 
Dr. Micheline C. Roufosse, Principal Investigator 


December, 1980 


During the reporting period ending on July 31, 1980, no work 
was performed on this project. Results of the work prior to this 
period were presented at the spring AGU meeting in Toronto, 
Canada, May 25-27. No project funds were spent during that 
period except in support of this presentation. 


During the last reporting period, from 1 August through 31 
October, most of the work that was done dealt with the northern 
section of the Mid-Atlantic Ridge and the Azores area, where we 
want to study the short and intermediate wavelengths in the geoid 
height spectrum. In order to accomplish this, we retrieved all 
the SEASAT and GEOS 3 satellite passes which crossed the area 
between latitudes 0 degrees and 40 degrees North and between 
longitudes 300 degrees and 360 degrees East. We retrieved 29 
SEASAT passes in that region but, unfortunately, many passes were 
almost repetitive; these 29 passes offered only 3 distinct 
groups of passes. Therefore, we decided to complement this data 
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Get with GEOS 3 data. The GEOS 3 data that were used have been 
corrected for bias and trend. The long wavelengths 04000 km) in 
the geoid heights have been removed by subtracting a reference 
geoid calculated with the GEM-7 coefficients up to degree and 
order 10. The geoid heights are in the process of being 
contoured using Mercator projections in collaboration with Drs. 
Parsons (Massachusetts Institute of Technology) and McKenzie 
(U. Cambridge, England) . By suppressing broad features of large 
magnitude, the process of filtering out the long wavelength 
features increases the resolution of the geoid obtained. These 
remaining geoid heights will then be studied in terms of the 
aging and cooling of the lithosphere. 

The small wavelength features in the geoid are studied on a 
track by track basis. The correlation between the geoid heights 
and the bathymetry is studied by using a thin elastic plate model 
with variations in crustal thickness. Typical crustal 
thicknesses along the North Atlantic Ridge which are derived by 
this method are of the order of 10 km, in agreement with 
observations in that region. 

Predictions are obtained in the Azores area using an Airy 
model and a thin elastic plate model. The results obtained so 
far tend to favor the Airy model with a crustal thickness of 
about 20 km. Further work is necessary, however, in order to 
gain more confidence in the validity of these preliminary 


results 


PAGE 3 


Estimated Fur ding Status 


Of the total amount of $50,000, approximately $29,600 
been expended. This represents about 60% of the fund; about 
of the total work has been completed. 


have 
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1 . 


In this work, we will study the behavior and mechanical 
properties o€ the lithosphere. This is a prerequisite to an 
understanding of the mechanisms and processes that occur in the 
earth's mantle, which are masked by the lithosphere. For this 
study, we will use the geoid heights derived from the GEOS 3 and 
SEASAT radar altimeters. The study of the correlation e;cisting 
between bathymetry and geoid heights gives information on the 
mechanical properties of the lithosphere such as its thickness, 
which is related to the age of the lithospheric plate. By 
probing several locations spanning varied temporal situations, we 
are able to retrace the time evolution of the lithospheric 
plates . 

We also believe that a systematic study of this type will 
lead to one of the following predictive capabilities; location 
of unknown seamounts from the knowledge of the geoid heights, 
prediction of gravity field from the knowledge of the topography 
and past history of the relief, or insight into the history of a 
region from the knowledge of both topography and geoid heights. 
To achieve this aim, several seamount chains, islands, and ridges 
must be investigated. Several geographical areas have already 
been studied, including the Walvis Ridge in the Southern Atlantic 
Ocean and several seamount chains in the Central Pacific Ocean. 
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2 . 

In this paper, we have restricted ourselves to the study of 
linear bathymetric features, i.e., seamount chains. In studying 
such features, it is legitimate to consider individual satellite 
passes crossing the chain of seamounts at an angle as close as 
possible to 90 degrees. It is therefore unnecessary to go 
through the tedious procedure cf adjusting all satellite passes 
into a coherent network by calculating and removing a bias and 
trend. Furthermore, dealing with individual satellite passes is 
ideal for studying the evolution along a seamount chain. 

To account for the geoid signals observed so far, two models 
were developed: the thin-elastic-plate model and the Airy model 
(see, for example, McKenzie and Bowin, 1976). T)ie 
thin-elastic-plate model applies when a load has been developed 
on a thin plate of finite thickness which is subsequently 
deformed. The magnitude and wavelength of tne deformed area 
depend mostly on the flexural rigidity, which is itself 
proportional to the cube of the lithospheric thickness. The 
Airy-type model applies when the load was created simultaneously 
with tae lithosphere (or on zero thickness lithosphere) and 
developed light roots in order to establish local or regional 
isostatic equilibrium. The methodology has been described in 
Roufosse (1979) and will be summarized here. 
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The first step in our procedure is to calculate theoretical 
filters in wavenumber space using variable parameters for both 
models (for example, see Figure 1) . We then Fourier-transformed 
these functions in normal space. Next, the bathymetry was 
reconstructed as rigorously as possible along the subsatellite 
position using available contour charts (the Uchupi charts were 
used in the Atlantic Ocean) . Lastly, the filters were convolved 
with the bathymetry in order to produce a sot of theoretical 
geoid profiles to be compared with the observed geoid heights. 


3. letoiii pgii.ofjD.gd Doling TJii£ Rep orti ng Egxlod 

Much of the work consisted of organizing and editing the 
SEAEAT radar altimeter data. The data are organized into 
individual files, each corresponding to a specific revolution 
number and further divided into sections corresponding to 
specific oceans. A program has been written to retrieve the 
satellite passes crossing a given geographical area. 

Hav/aii, in the Pacific Ocean, was chosen as the test area. 
An extensive study of this area, using the GEOS 3 radar altimeter 
data, has been done at SAO in collaboration with Dr. B. Parsons 
(Roufosse, M.C., Parsons, B.E., EOS Trans., A.G.U., Jlfi, 1213, 
1977) . Ninety SEASAT passes have been retrieved, plotted, and 
edited in this region in order to assess, in a known case, the 
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gains in accuracy offered by SEASAT, and therefore indicate which 
steps should be taken next. GEOS 3 has been proven to be 
inadequate for studying features that offer a weak geoid signal, 
such as the East Pacific Rise, because the noise level was too 
high. Figure 1 shows a comparison between the geoid heights 
derived from GEOS 3 and SEASAT over the Hawaiian Seamount chain 
after the SEASAT data has been edited. The obvious gain in 
accuracy obtained and the low noise level observed with SEASAT 
will make it possible to study such interesting features as the 
East Pacific Rise as soon as the Southern Hemisphere data have 
been received. 

The SEASAT data received did not prove to be fully adequate 
for the work proposed; therefore, we used a combination of 
SEASAT and GEOS 3 data. 

The two areas that look most promising in the SEASAT data 
set available at this time are the New England seamounts and the 
northern section of the Mid-Atlantic ridge. 


Two 

typic 

al passes over 

the 

New England seamounts 

can 

be 

seen on 

the 

top profiles 

of 

Figures 2 and 3. The 

signals 

observed 

have 

a wavelength of 

150 

to 200 km and an intensi 

ty 

of 2 

to 4 

m. 

It is similar 

to 

the signal observed 

in 

the 

Marshall- 

-Gil- 

.'t seamounts area. 
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Tho correlation study between the bathymetry and the gcoid 
height yields information on the mechanical properties of the 
lithosphere; in this case, all the values found for the flexural 
rigidity are very low. These values range between 1 x 10^® and 4 
X 10^^ dyne-cni. The lithosphere is known to be very old in this 
region and therefore it is very strong. In order to account fo*- 
the low values of the flexural rigidity, wo had to postulate that 
the seamounts were very old and, therefore, were formed where the 
lithosphere was younger and thinner. The information gathered on 
the ago of the seamounts confirms that hypothesis. 

We also studied the Rio Grande Rise in the South Atlantic 
Ocean. This area was chosen in order to be compared with the 
previously studied Walvis Ridge area and because it is located 
symmetrically with respect to the Mid-Atlantic Ridge. The geoid 
signals present in the Rio Grande Rise region are weak and broad 
and mainly resemble those observed in the eastern and central 
sections of the Walvis Ridge. At this point in the study, we 
believe that the Rio Grande Rise has been formed simultaneously 
with the eastern and central sections of the Walvis Ridge by a 
hot spot located on the Mid-Atlantic Ridge. An Airy-type model 
of compensation best describes the signal observed. Because the 
SEASAT data in the Southern Hemisphere was not available until 
very recently, the study had to be complete'; with GEOS 3 data for 
good data coverage. 
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Now, we are in the process of retrieving several passes over 
the northern Mid-Atlantic Ridge and over the Azores Islands in 
order to study both the short and intermediate wavelengths that 
are present in the geoid spectrum. 

Twenty-nine SEASAT pasyes and sixty GEOS 3 passes have been 
retrieved in that area. The GEOS 3 data are being filtered for 
long wavelengths (a reference geoid of degree and order 10 
calculated with GEM 7 is removed from the raw data ) , and 
corrected for bias and trend using Dr. Rapp's information. It 
will then be plotted on Mercator projection in a collaborative 
effort • with Dr. Parsons (MIT) and Dr. McKenzie (Cambridge 
University, U.K.). This map will give information on the 
intermediate wavelengths contained in the geoid and will be 
suitable for a time dependence study. The SEASAT data are 
studied on a track-by-track basis, and attention is focussed on 
the small wavelength features. We are presently in the process 
of retrieving the bathymetry along individual passes using all 
available information. 
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Figure 2. The top profile represents the observed geoid heights for the 
Seasat revolution number 529 over the New England seamounts. 

The bottom profile represents the predicted geoid heights 
calculated with a value for the flexural rigidity of 1.3 x 10^- 
dyne-cm. 







